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I . Introduction 

This report describes a two dimensional stochastic process 
and a program for its simulation. The process can be considered 
to describe a random track of predetermined length that origi- 
nates and terminates at predetermined end points. The process 
is a modification of a px'ocess that is without the length and 
termination constraints. This process, which is also described, 
is a generalization of the random tour process described by 
Washburn in Reference 1. 

A particular process is defined by the track length, the 
maximum expected track segment length and the distance between 
the end points. For motion at a constant speed, the track 
length constraint is equivalent to a transit time constraint. 

In Section II, the process and the unmodified process are 
described and some of their characteristics are discussed. In 
Section III, the simulation program is described. By using 
the program, graphs and daoa can be produced that represent 
realizations of the process. The program, which is written in 
BASIC, is for an HP“85. 
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II. A Description of the Process 

The process, as it is described here, generates a random 
track of predetermined length that connects two predetermined 
end points. The track originates at a predetermined starting 
point, the first point, and terminates at a predetermined stopping 
point, the final point. The track is continuous between these 
two points and consists of a sequence of straight line track 
segments. Some representative tracks are shown in Figure 1. 

A sequence of rectangular coordinate systems is used to 
define the process. The first is a reference system for the 
first track segment, the second is a reference system for the 
second track segment and so on. The x-axis of the first system 
is coincident with the line joining the first and final points, 
the origin is at the mid-point of the line and the system is 
oriented so that the final point is on its positive x-axis. 

For each successive system, the x-axis is coincident with the 
line joining the last point of the preceding track segment and 
the final point, the origin is at the mid-point of the line 
and the system is oriented so that the final point is on its 
positive x-axis. Each track segment is determined by two quan- 
tities: an angle 9 and a length r. The angle 9 is measured 

clockwise from the direction of the positive y-axis of the 
reference coordinate system of the track segment. It is the 
relative course that will generate the track segment. The 
length r determines the track segment length. 
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Figure 1. 



Some representative tracks. The tracks were generated 
by using the program that is described in Section III. 
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The process is defined as follows: For each track segment 

except the last a relative course 9 is chosen from a distribu- 
tion with density function 

^ 9 ( 9 ) “ (1/2tt) b/ (a - c sin 9) 

where 0 < 9 < 2 tt, 2c is the distance between the first and 
final points, 2a is the track length with a > c > 0 and 
b = (a - c ) ^. Next, a length r is chosen from a conditional 
exponential distribution with density function 

f^|g(r|9) = [ 1 / 6 (9) ]exp[- r/6(9)] 

where 0 < r, 5(9) = (a - c)6/(a - c sin 9) and 5 is the maximum 

expected track segment length. Using the parameters a^,, b^^ and 

2 

Cj^ defined below, if r > 2a^^, then r^^ = b^ /(aj^-c^ sin 9j^) . In 
this case, the (k+l)st track segment is the line that joins the 
last point of the kth track segment to the final point and the 
process terminates. If r < 2a^, then r^ = r. In this case, the 
last point of the kth track segment and the final point deter- 
mine the (k+l)st reference coordinate system and the process 
continues. The parameters in the above expressions are defined 
as follows: 2aj^ is the remaining track length with 2a^^ = 2a 

and 2a, = 2a, , -r, , for k > 1; 2c, is the distance from the 

k k-1 k-1 k 

last point of the (k-l)st track segment to the final point with 
2c^ = 2c and b^^ = (a^^ -c^^ ) The program that is described in 
Section III can be used to simulate this process. 

Figure 2 illustrates the generation of a three segment 
track. The ellipse on which the second segment terminates 



4 





Figure 3. The first rectangular coordinate system, the 
associated polar coordinate system and the 
reference ellipse. 
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and the third segment begins has one focus at the last point of 

the first segment and the other at the final point. With 

the distance between these two points, 2 a 2 = 2a - r^ the 

2 2 2 

remaining track length and b 2 ~ ^2 ~ ^2 ' ellipse is 

2 2 2 2 

defined through x /a 2 + y /b 2 = 1 in the second reference 

2 

system. It is also defined by r = b 2 /(^2 ~ ^2 

polar coordinate system with its origin at the last point of 

the first segment and its polar axis parallel to the y-axis of 

the second reference system. 

As a basis for discussing the relative course and track 

length distributions that are associated with the process, 

consider the polar coordinate system with its origin at the 

first point and its polar axis parallel to the y-axis of the 

first reference system. Next, consider a reference ellipse 

with one focus at the first point and the other at the final 

2 

point that is defined by r = b /(a - c sin 9) . The coordinate 
system and the ellipse are illustrated in Figure 3 . 

The distribution of the relative course 9 and the condi- 
tional exponential distribution of the length R are both deter- 
mined by p(9) = (a-c)/(a - c sin 9), the ratio of the value of 
the radial coordinate of the ellipse at 9 to its maximum value 
at 7 / 2 . Note that when a track segment's relative course is 
tt/ 2 the segment is directed toward the final point. The dis- 
tribution of the relative course G is determined by the require- 
ment that the ratio of the probability that a relative course 
is between 9 and 9 + A9 to A9 is proportional to o(9) . The 
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conditional exponential distribution of the length R is deter- 
mined by the requirement that the ratio of o(9), the expected 
track segment length for 0, to 6, the expected track segment 
length for tt/2, equals p(9). 

If the first point and final point are coincident, then 
c = 0 and the reference ellipse becomes a circle. In this case, 
the distribution of the relative course 0 is uniform with den- 
sity function ^0(9) = 1/2tt where 0 < 0 < 2it, the relative course 
0 and the length R are independent and the distribution of 
R is exponential with density function f„(r) = (1/5) exp {-r/5) 
where 0 < r. 

Consider a process that is defined in the same way except 
that the length of a track segment is always chosen to be equal 
to the value determined for R. For this process, a is not the 
track length, there is not a length constraint, and the final 
point is not an end point, there is not a termination constraint 
The final point is a reference point of the process, literally, 
focal point of the process. If c = 0, the process is identical 
to the random tour process described by Washburn in Reference 1. 
If c > 0, it is a generalization of that process. 

By inspection of the density function, it can be seen that 
both the mean and the mode of the relarive course i are equal 
to ’T/2. For the constrained process, this implies that the 
most likely next to last and last track segment relative course 
values are tt/2 and -"t/2. In Appendix 1, it is shown that the 
expected value of the length R is ao/Ca+c) . For the uncon- 
strained process, this implies that the average track segment 
length is ao/{a-^c). 



A new process could be generated by redefining the para- 
meters a and c. For example, one could use the definition 
a = 2c where c is again the distance between the first and final 
points. In this case the degree of focusing of the process 
would be constant. However, relative to the process described 
here, it is conjectured that, for track lengths greater than 2a, 
the process would concentrate track points about the final point, 
for track lengths less than 2a, the process would concentrate 
track points about the first point. In the simulation program, 
for a track of m segments where m > 2, the maximum length of 
the first m-2 track segments is computed. And, it is also con- 
jectured that, for a track length less than 2a, the probability 
that the length of the last track segment would exceed this 
statistic would be greater than that for the process described 
here. 
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III. The Simulation Program 

The simulation program requires three inputs: The distance 

D between the end points, the track length TL and the maximum 
expected track segment length 6. The program is for an HP-85 
and is written in BASIC for that computer. The program listing 
is in Appendix 2. After loading the program and initiating RUN, 
a user is prompted for these inputs and advised of their limits. 
An example of the prompting with a set of inputs is shown below. 



D I S T M H C E 0 i.1 0 ? 

1 0 

TRRCK LEHGTH D <TLi.50 ;'*• 

50 

0 E L T ,R 1 . '> ’? 

1 

i. 

Table 1. The input limits are based on display and memory 
requirements . 

As a simulation progresses, the developing track is displayed 
on the HP-85 CRT. After a simulation is completed, a user is 
given the option of displaying track data or printing the track 
and the track data. In the CRT display of a -rack, uhe starting 
point is on the left side and the stopping point is on the right 
side. In the printed output, a track is rotated 90° so that the 
starting point is at the top. Figure 4 shows sections of the 
printed output. In the columned data, I is the track segment 
number, f is the track segment course for a starting point that 
is due west of the stopping point, R is the track segment length. 
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0 is the track segment relative course and X and Y are the 
coordinates of the last point of the track segment in the first 
coordinate system. If the starting point is considered to be 
due west of the stopping point, then the positive x-axis of the 

I 

first coordinate system is directed east and the positive y-axis 
is directed north. Recall the origin of the first coordinate systenij 
is at the mid-point of the line joining the starting and stopping j 
points. ! 

For the track shown in Figure 4, the termination criterion 
was satisfied at the simulation of the 49th track segment. In 
the first line below the columned data, 0 is the average track 
segment relative course where -90 °< 9 < 270°, R is the average 
and RM is the maximum track segment length for the first 49 
track segments. In the second line, <I>(51) is the relative course 
of the 51st and last track segment and R(51) is its length. 

If a track simulation reaches 97 segments without termina- 
tion and the termination criterion is not satisfied for the 98th 
segment, the simulation is stopped. In this case, the values 
on the first line below the columned data refer to all of the 
98 track segments and "TRUNCATED" is on the second line. For 
both cases, the average relative course 9, the average track ' 

length R and the maximum track length RM relate to the uncon- 
strained stochastic process. j 

Figures 5, 6 and 7 show sections of the printed output for ^ 

three additional simulated tracks. The plots illustrate the 
effect of the value of 5 on the plot scale. Doubling the value 
of 0 doubles the distance between points on the plot. This 
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increases the likelihood that the HP-85 CRT will display a com- 
plete track. The plots illustrate, to some extent, the increased 
focusing that occurs as the value of TL approaches the value of 

D. 

A second simulation program to generate track data in a non- 
graphical but more efficient manner was written. A preliminary 
analysis of some data generated by the simulation has been made. 

/N 

In Table 2 below, some estimates p with associated 95% confidence 
intervals are given for p, the probability that the last track 
segment length will exceed the maximum track segment length of 
the preceding track segments. 

D TL 6 p 95% confidence 

interval 



10 


50 


1 


.60 


.57 


< 


P 


< 


. 63 


10 


50 


2 


.55 


.52 


< 


p 




. 5 8 


10 


50 


5 


.53 


.50 


< 


P 


< 


. 57 


10 


20 


1 


.44 


.40 


< 


P 


< 


.47 


10 


20 


2 


.43 


. 39 


< 


p 


< 


.46 


10 


20 


5 


45 


.42 




P 


< 


.48 



Table 2. Estimates p for p, the probability that the last 

track segment length will exceed the maximum track 
segment length of the preceding track segments. 
Each estimate is based on a sample size of 900. 
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. A program generated rrack. The track dat.a for 
segments 4 through 47 have been omitted. 



12 



Figure 4 
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DELTR= 1.00 




Figure 5. 



A program generated 
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igure 6 . 



A program generated track. 
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D I S T R H C E = 1 0 0 ij 



Figure 7. 



TRRCK LEHGTH= 20 . 00 
DELTR= 2 00 




A program generated track. 



Appendix 1 Some Mathematical Relationships 



In the simulation, values of 0 are generated by using the 
inverse transform method. This method requires the use of 
the cumulative distribution function Fg(9) and it can be found 
as follows: For -tt/2 £ 9 < tt/2 , 

9 

F„(9) = / (1/2tt) [b/(a-c sin 9)] d9 

® -tt/2 

= 2 [b/ (a+c) ] { tan (a tan (9/2)-c)/b] - tan ^[-(a+c)/b]} 

For tt/2 < d < 3n/2, Fq(9) = 1 - Fg(7T-9). 

The expected value of R for the unmodified process that is 

given in Section II can be found as follows: 

tt/2 

E(R) = 2 / E(R]0=9) fg(9) d9 

-r/2 

-/2 

= / iS/^) b (a-c) [l/ (a-c sin9)^] d9 

-tt/2 

= ( 26/tt ) [a/ (a+c) ] { tan ^[b/(a+c)] + tan ^[(a+c)/b]} 

= ao/a+c 

If the termination criterion is satisfied on the first 

track segment, the track will consist of only two track segments. 
2 

With r^(9) = b /{a-c sin 9), the probability P 2 of this event 
can be found as follows: 
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tt/2 

2 / P(R > r^(9) |0 = 9) fg(9) d9 

-tt/2 

tt/2 

2 / exp[-r^(9)/o (9) ] f^iQ) d9 

-tt/2 



exp[- (a+c) /5] 
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The HP-85 Simulation Program 



L ’•I' C 0 
0 y > 



1 0y ) R ( 1 00 . ' ■ T ' • . 1 00 
i 0 0 f 



\ 1 



20 D I Sr " D I S'^RNCE “ 8 '' CHP $ 

S S > “ D " ; C H R $ < 1 S S > ; " 1 8 > " ; 

30 INPL'T OS Cl=C'-'2 
40 IF O>10 OR D<0 THEN 20 
50 OISP “TRRCK LENGTH".." 

RJ < 60 > ; “ TL •' CHR^ •: 1 SO ; 

60 INPUT L 
70 IF L.>50 OR L-:=0 THEN 50 EL 
Rl=L.--2 

S0 D I SP " OELTPi " .1 " 1 " CHR-j: t 1 S 

CHR S < 1S> : " > ” 

20 INPUT Di 
10101 I.F Ol< = 0 THEN 90 
110 RRNOOmZE 

120 h 0i> =-C 1 S Y<0>=0 S F":. 0i.>=0i 
fi; ( 0 .1 = 0 6 T <. 0 > = 0 
130 81=0 S B2=0 S B3=0i S Rl=0i 
fi=Hl S C = C1 S S=0i S V = 0i S 
9 

1R0 GCLERR S E=10 t:SQR(;Oi ; 

1 50 SC RLE -E E .. - . 75;icE :• . 75:? E 
160 MOVE C1..0’ e L.R8EL "o" 

170 MOVE -Ol.. 01 S LRBEL "o" 

190 FOR 1=1 TO gP 
120 U=RHD S G1=0 
2 01 0 IF U •; = . 5 THE f i 2 2 0 
210 U=U- . 5 S Gl=l 
220 HI =SQR t; < Hi 1 > R 1 -C I S 

= S Q R ( 1 - C 1 1 C 1 ••• ■ (. R 1 1 R 1 ' 

230 T = 2-«RTN < H2 '^:TRN < P I .t:iJ-RTN N 1 
+ C 1 •••'R 1 ? 

24y Gl=l THEN T=PI-T 
2 5 0 P = <. R 1 - C 1 ■' . ‘ R 1 - C 1 1 S I f < ( I > > 
260 D2=D1 :rP 
2 70 R = - D 2 L 0 G C 1 - R N 0 > ) 

290, F = T H 

2 ? 01 G 0 S U B 9 0 0 

300 i~ h2>C2 then 340 

3 1 0< R = ( M -t R - ;« C > .•••' '■ R - C t S I H T > 
3201 GOSUB 80101 

3 4 01 >G R 4 ‘:5) I* T *;i V r' ^ i 

3 5 01 R = R 2 S C = C 2 S 8 = - R T N 2 Y .■ X > 

C* ( T j ix' 

3 i H = T G 0 S L‘ B ”” 7 k? 

3 ? 9 T »■’ I > = H 

3B9 H = F GOSUB 779 

3 ’? r I > = H 

*^99 i i“ E* 2 “ 1 T H £ H *4 3 

>^19 R1 > 

4 29 W = ’,.j T — T . I l»J = l«4 -f- I' Jr. — W > ■ I 

439 ne:<t I 

440 o;;p \i DI3P '* "RUHChTED ‘‘ 2 

I 9 E 3 = 1 2 G 0 S i J B S' 5 9 2 G 0 T 

4 > VI 

4 ^ VI Q ►»■' v*4 _ 1 •• 

4 S’ 9 K =1-^1 2 = 2 C 2 7 f K > = 1? 9 

H = 3 P I ■ 2 2 iG 0 S- U b 779 2 ^ . 
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520 
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560 
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570 
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530 
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650 
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END 
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